arXiv: 1504.0557Ivl [math-ph] 21 Apr 2015 


Vector Riemann-Hilbert problem with almost periodic and 
meromorphic coefficients and applications 


By Y.A. Antipov 

Department of Mathematics, Louisiana State University 
Baton Rouge LA 70803, USA 


Abstract 

The vector Riemann-Hilbert problem is analyzed when the entries of its matrix 
coefficient are meromorphic and almost periodic functions. Three cases for the mero¬ 
morphic functions, when they have (i) a finite number of poles and zeros (rational 
functions), (ii) periodic poles and zeros, and (iii) an infinite number of non periodic 
zeros and poles, are considered. The first case is illustrated by the heat equation for 
a composite rod with a finite number of discontinuities and a system of convolution 
equations; both problems are solved explicitly. In the second case, a Wiener-Hopf 
factorization is found in terms of the hypergeometric functions, and the exact solu¬ 
tion of a mixed boundary value problem for the Laplace equation in wedge is derived. 
In the last case, the Riemann-Hilbert problem reduces to an infinite system of linear 
algebraic equations with the exponential rate of convergence. As an example, the 
Neumann boundary value problem for the Helmholtz equation in a strip with a slit 
is analyzed. 


1 Introduction 


Many physical models described by boundary value problems for elliptic, hyperbolic, and 
parabolic equations reduce to a system of n convolution equations on a finite segment. 
Such systems are equivalent to a vector Riemann-Hilbert problem (RHP) with a block- 
triangular matrix coefficient 


G{a) 


e*"I 0 \ 

g{a) e-*“/ )’ «e(-oo,+oo). 


( 1 . 1 ) 


where I is the order n unit matrix and g{a) is an re x re matrix. The functions e*“ 
and are almost periodic functions (Levitan, 1953), and their indices are infinite: 
inde*" = -boo and inde“*“ = —oo. If re = 1 and g{a) is a rational function, then the 
RHP admits a closed-form solution (Ganin, 1963). In the general case of the function 
g{a), even when re = 1, there is no procedure for solving the RHP with the coefficient 
(ll.ljl in closed form. Novokshenov (1980) analyzed the singular convolution equation 
/g“[(x — t)~^ + k{x — t)]u{t)dt = f{x), 0 < X < a, k{x) € L 2 {—oo, -boo), and showed that 
its solvability and the representation formulas for the solution are obtained in terms of 
the Wiener-Hopf factors of the associated matrix (ini) (re = 1). A theory of factorization 
of matrices (II.Ih when g{a) is an almost periodic function (including an almost periodic 
polynomial) was developed by Karlovich and Spitkovsky (1983), Spitkovsky (1989) and 
their coauthors (see for example Bottcher et al, 2002). 

In many applications to physical models, g{a) is a meromorphic function. Antipov 
(1987, 1989) considered contact problems on an annular stamp and reduced it to a RHP 
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with the matrix coefficient p.ll) . g{a) = ^r(^)r (7 +1 — f )[r(| + ^)r (7 +1 — f )]“^, 7 = 
0,1. The problem was transformed into an infinite system of linear algebraic equations of 
the second kind with the exponential rate of convergence and solved in terms of recurrent 
relations. Recently, this technique for employed for the solution of an integro-differential 
convolution equation arising in fracture with surface effects (Antipov and Schiavone, 
2011 ). 

Another class of matrices of the form 


G{a) 


gii{a) e*"fi(i 2 (a) \ 
e"*" 5 ' 2 i(a) 522(a) j ’ 


a E (—00, +00), 


( 1 . 2 ) 


needs to be factorized to solve the systems of convolution equations kji * (pi + kj 2 * 4>2 = 
fj{x), ttj < X < 00 , supp^j C [aj,oo), j = 1,2, ai = 0, a 2 = a. A method of inte¬ 
gral equations for factorizing matrices of such a structure was proposed by Abrahams 
& Wickham (1990). Onishchuk (1988), Antipov k. Arutyunyan (1992), Antipov (1995, 
2000) worked out a technique for vector RHPs when 57 (a) are meromorphic functions, 
and applied it to static fracture and contact problems. This method ultimately requires 
solving inhnite systems of linear algebraic equations of the second kind with the expo¬ 
nential rate of convergence. 

In this paper we aim to develop further the methodology for the RHPs whose matrix 
coefficient entries are meromorphic and almost-periodic functions and apply it for the 
solution of some model physical problems. In Section 2, we consider the heat equation 
for an inhnite rod ut = a?‘{x)uxx + g{x,t) with a piece-wise constant diffusivity a{x). 
Deconinck et al (2014) applied the Fokas method (2008) to hnd an exact solution in the 
homogeneous case, (7 = 0 , when the diffusivity is a piece-wise constant function, and 
(i) the rod is hnite and the diffusivity has one or two points of discontinuity, and (ii) 
the rod is inhnite and the function a{x) is discontinuous either at 0 and 00 , or at two 
hnite points and inhnity. In particular, for an inhnite rod with a{x) = a±, > 0, 

and (7 = 0 , they managed to represent the solution by double integrals with one-sided 
Fourier internal integrals and contour external integrals over the boundary of inhnite 
wedge-like domains. In our case, the rod is inhnite, and the diffusivity has n -|- 1 points 
of discontinuity. We show how the problem can be reduced to an order-n RHP and 
solved exactly in terms of the solution of an associated hnite system of linear algebraic 
equations. In particular, when n = 1, and the points of discontinuity are 0 and 00 , 
we simplify the solution and derive an exact formula for the temperature. If g( = 0, 
then the solution is given by a sum of two one-dimensional integrals over the intervals 
(— 00 , 0 ) and (0, 00 ). The formula found generalizes the classical Poisson formula for the 
homogeneous inhnite rod. We also briehy describe another approach for the solution of 
the heat problem that employs the Laplace transform and the theory of discontinuous 
one-dimensional boundary value problems. The solution constructed by this alternative 
approach coincides with the one found by the RHP technique and is applicable to both 
cases, hnite and inhnite, and for any hnite number of discontinuities. At the end of 
Section 2, we present an alternative approach for solving the Abrahams-Wickham system 
of convolution equations (Abrahams k Wickham, 1990) by reducing it to the case (11.21) 
with gij being rational functions. Notice that their procedure for the RHP hinges on a 
solution of an auxiliary system of integral equations, while our approach bypasses extra 
integral equations. 

In Section 3, we propose a factorization method for matrices (ll.2jl when the functions 
gij are meromorphic and their zeros and poles are periodic. The method is illustrated 
by solving a mixed boundary value problem for the Laplace equation in a wedge. This 
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problem for a particular choice of the boundary data was analyzed by Abrahams & 
Wickham (1990). They rewrote the problem as a vector RHP with the coefficient of the 
form p.2l) and reduced it to a system of two integral equations to be solved numerically 
by an approximate method. The approach we propose does not require solving any 
auxiliary integral equations. It derives Wiener-Hopf factors of the matrix coefficient in 
terms of the hypergeometric functions and ultimately yields a closed-form solution of the 
physical problem. 

In Section 4, we generalize the method for RHPs with the matrix coefficient of the 
form (II.Ih for the dynamic case. As an illustrative example, we take the Nemann bound¬ 
ary value problem for the Helmholtz equation in a strip with a finite slit. This problem 
describes antiplane strain deformation of a strip when the strip boundary is free of trac¬ 
tion, and the Mode-HI crack faces are subjected to oscillating loading. At the same 
time it could be interpreted as a model of sound transmission in a waveguide when an 
acoustically hard finite screen is placed inside the waveguide. We derive the solution by 
quadratures and some exponentially convergent series with the coefficients determined 
from a rapidly convergent infinite system of linear equations of the second kind. 

2 Matrices with almost periodic and rational entries 

In this section we derive an order-n vector RHP associated with the heat equation in a 
rod with a piece-wise constant diffusivity and conductivity and show that it admits a 
closed-form solution. To verify the procedure, we derive the solution by the standard 
technique of discontinuous one-dimensional boundary value problems that bypasses the 
RHP. To give an extra example, we solve a system of two convolution equations. 

2.1 Heat equation with piece-wise constant coefficients: RHP for an 
order-n vector-fnnction 

The problem under consideration is one of heat conduction for an infinite rod with a 
piece-wise constant diffusivity o^(x) = k{x)/\cpp{x)] 

ut = a^{x)uxx + 9 {x.,t), |x| < oo, x 7 ^ 6o, 6i,..., 6n-i, t > 0. 

^\x=b~ ~ ^\x=b+^ ^j^x\x=b~ ~ ^3+^^^\x=b+' J = 0, 1, . . . , n — 1, t > 0, 

3 3 3 3 

u\t=Q = /(x), |x| < oo, (2.1) 

where u{x,t) is the temperature, g{x,t) = {cpp)~^go{x,t), go{x,t) is the heat source 
density, /(x) is an initial temperature, k{x) is the thermal conductivity, p{x) is the 
density, Cp is the specific heat capacity, and 

a(x) = aj > 0, k{x) = kj, p{x) = pj, x € (6j_i,5j), j = 0,1,... ,n. (2.2) 

Here, we assumed 6_i = —oo and bn = -|-oo. In what follows we reduce this physical 
problem to an order-n vector RHP. Introduce first the Laplace transforms 

PCX) PCX) 

u{x]p)= / ue~^^dt, g{p) = / ge~^^dt, Rep>0, (2.3) 

Jo Jo 

and obtain from (12.ip a discontinuous one-dimensional boundary value problem. It reads 
0?{x)Uxx - pu = -f - g, |x| < oo, X / 
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'^\x=b~ ~ '^\x=b+’ ^j'^^\x=b~ ~ ^j+^'^^\x=b+^ j — Oj 1) • • • j''T- 1- (2-4) 

3 3 3 3 

To apply further the two-sided Laplace transform, we introduce new functions of the 
parameter p 

and integrate by parts 

n—1 


Pojip) = 

j3ij{p) = UjUxl^^f,- - a]+^Ux\^^f,+ , j = 0,1,... ,n - 1, 


/ OD ^ POO 

{x)uxxe~^^dx = ^(/3ij + s/3oj)e~^^^ / a^{x)ue~^^dx. 

-OO j_Q J — oo 

We split now the integral in (|2.6I) into n -|- 1 parts and denote 

rbo 


( u{x]p)e ^^dx = e ^^°Uq{s), 

J —OO 


(2.5) 


( 2 . 6 ) 


Oj-l 


where 


u{x\p)e ^^dx = e (s) = e (s), j = 1, 2,..., n — 1 

POO 

/ u{x-,p)e~^^dx = {s), 


/ O poo 

u{x + bo]p)e~'''^dx, Uq{s)= / u(x + 6n_i;p)e“®^dx, 

-OO Jo 


(2.7) 


Utis) = 


0 

[ u{x+ bj;p)e~^^dx, U~{s)= [ u{x + bj-i]p)e~^^dx, 

Ibj—i—bj Jo 


j = ( 2 . 8 ) 

The functions {s) and U^{s) {j = 0,1,..., n — 1) are analytic in the domains D~^ = 
{Re s < 0} and D = {Re s ^ respectively. We emphasize that except for Uq (s) all 
the other functions Up{s) (j = 1 , 2 ,... , n — 1 ) are entire functions in the half-planes D^. 
In these notations, the one-dimensional boundary value problem (j2.4p can be recast as 
the following order-n vector RHP: 

UUs) = j = 1, 2 ,... , n - 1 , 


71—1 


'^mj{s)e {s) + mn{s)e (s) = H{s), s £ L, (2.9) 

j=o 

where L is the positively oriented imaginary axis (the domain D'^ is on the left), mj{s) = 
a^jS^ — p, j = 0,1,... , n, and H{s) is given by 


POO ^ 1 

H{s) = - [f{x)+g{x-,p)]e~"^dx - ^ s/3oj). 

j=0 


( 2 . 10 ) 


The vector RHP (12.9p can be transformed into a finite system of linear algebraic equa¬ 
tions. To do this, without loss of generality, we assume 6o = 0 and rearrange the RHP 
(1^ as follows: 


Uois) + ^U^{s) + ... + + ... + e-^^"-^t/-_i(.) 


mo 


mo 


mo 
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mo mo 


m,- m, m,' 


m.- nriA 


n — 2, 


e"^"-i-^C/o+(s) + ... + + ... + e^(bn-i-6n-2)!I^[/+ (g) 

rr 7 „ 1 u \ J ^ _, n z\ / 


rrin-i 


rrin-i 

rT4- / \ TT-/ \ e^^"-^H(s) 

+ C^„_i(s) H - Uq (s) = -, s £ L. 

rUn-l rUn-l 


rrin-i 


( 2 . 11 ) 


Notice that the functions (s) {I = 0,1,... , j — 1) are analytic in and decay 

exponentially as s —>• 00 in D'^, while the functions , I = j + 1,... ,n — 1 

{U~{s) = Uq {s)) are analytic in D~ and decay exponentially as s ^ 00 , s G D~. On 
factorizing the functions 


mj+i{s) _ (s) 

mj(s) K-{s)’ 


j = 0,1,... ,n- 1, 


where 

K-{s)= + 

“i'S - VP aj+is + Vp 

and substituting this into equations (j2.11jl we have 



Re y/p > 0, 


( 2 . 12 ) 


(2.13) 


Utis) 


+ 


K^is) ’ K;{s) 

^7+1 (") 1 

K-{s) KHs) 


e^^^—U^{s) + ... + 

m,- nii 




sibj-bj+i)!!h±lu-{s) + ... + —[/o-(s) 

rrii rrii ^ 


-njis) 


j = 0 , l,...,n-l. 


(2.14) 


Here, U-^is) and T-L~{s) provide the splitting of the functions H{s)[mj{s)K^{s)\~^ 
into analytic parts in the domains D'^ and D~, respectively. In general, they are defined 
by the Sokhotski-Plemelj formulas 


n^{s) = ± 




e^^iH{a)da 


+ 


2mj{s)K^(s) 2TTi Jl mj{a)K^ (o')(o' — s) ’ 


s G L>, 


(2.15) 


and the Cauchy integral is explicitly evaluated by the theory of residues. Alternatively, 
this splitting can be obtained by representing the function H (s) as a sum of n integrals 
similar to (j2.8|] and then removing the poles. The second approach will be employed in 
the scalar case in section ESI 

Now, since the poles of the left- and right-hand sides are known, we apply the Liouville 
theorem. Crucial to the success of the method is the fact that the known functions in 
the system (I2.14h are meromorphic functions having a finite number of poles. That 
is why the left- and right-hand sides are rational functions with prescribed poles and 
unknown coefficients. These coefficients can be fixed by the requirement that the final 
solution Up{s) of the vector RHP has to have removable singular points at the poles 
lying in the half-planes D^. These conditions form a finite system of linear algebraic 
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equations for the unknown coefficients. The case of an order-2 RHP when the functions 
are meromorphic and have an infinite number of periodic and not periodic poles will be 
considered in sections [3] and 01 respectively. The functions /3oj{p) and f3oj{p) introduced 
in ()2.5I) are fixes by the conditions (|2.4I) . 


2.2 Generalization of the Poisson formnla for an infinite piece-wise 
homogeneons rod 

To clarify the procedure of finding the functions f3oj (p) and /3ij {p) , we consider the heat 
equation for an infinite rod composed of two semi-infinite rods having different constant 
diffusivities and conductivities. We assume that the initial temperature, fo{x), does not 
vanish at ±oo that is 


fo{x) ='y-9{-x)+ -f+9{x) + f{x), /(x) G Li(-oo,oo), (2.16) 

where 7 ± are nonzero constants and 9{x) = 1, x > 0 and vanishes otherwise. It is 
convenient to split the temperature, uo{x,t) (|x| < oo, t > 0), as 

uo(x, t) = u{x, t) + 7 _ 6 '(—x) -I- 7+6'( x ) (2-17) 

and determine u{x,t) as the solution of the boundary-value problem 

ut = a‘^{x)uxx + 9 {x,t), |x| < oo, x 7 ^ 0 , t > 0, 


u\x=Q— ^|a:=0+ T) k—Ux\x=0— |x=0+) ^ ^ Oj 

u\t=o = f{x), |x| < OO, (2.18) 

where a(x) = a_, x < 0, a(x) = a+, x > 0, and 7 = 7 + — 7 _. The problem is equivalent 
to the following scalar RHP: 




— p 


s‘^a?_ — p 


(2.19) 


where 

n dP o 

Pj = a^_^u{0~;p) - a+^u(0+;p), j = 0,1, 

POO pO 

= / e~^^udx, U^{s) = / e~^^udx, 

Jo J—CO 

rco pO 

H~{s)= e~^^[f{x) + g{x;p)]dx, H~^{s)= e~''^[f{x) + g{x-,p)]dx. (2.20) 

JO J —CO 

The functions and are analytic in the half-planes . The coefficient of the RHP 
is a rational function, and the functions U^{s) and U~{s) are recovered in the standard 
manner. 


U+{s) = 


sa- - y/p 


a^h- 


y/pia+ + a-) 


— An¬ 


sa- + 




U-[s) 


1 

sa+ + y/p 


a_/i+ 

y/p{a++a-) 


Ai- 


1 


sa+ - Vp 



a+h-{sa- + y/p) 
y/p{a+ + a-) 


( 2 . 21 ) 
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where 




Ao = 


a_ 

Po^ + a+Pi 
y/p(o+ + a_)’ 




= 


a+. 

Poy/P - Q-A 

y/p{a+ + o_) ’ 


( 2 . 22 ) 


Notice that the points s = G are removable singularities of the functions 

U^{s). The derivation of representations for the functions (3o{p) and f3i{p) requires 
inversion of the Laplace transforms in (12.210 . This implies 


+ 


u{x;p) = 
1 


2a_/i+ — (a+ — a-)h- Ai 


p-\/pa;/a+ 


‘^a+y/P Jo 

u{x‘,p) = 


2a+{a++ a-)y/p a+ 

roo 

/ + 0 <x<oo, 

Jo 

g^x/a_ 


2 a+/i_ + (a+ — a_)/i+ ^ 


2a_(a++ a_)^/p a_ 

+ /_J/(6 -oo < X < 0. 


(2.23) 


The function u{x;p) and its derivative Ux{x]p) are discontinuous at the point x = 0 and 
due to (j2.18jl have to meet the conditions 

u\x=o~ ~ ^U=o+ ~ k—Ux\x=o~ ~ k-\-Ux\x=o+ ~ 0) (2.24) 

On satisfying these conditions we eventually determine the functions /3o and /3i as 
1 


/3o = 


o+A;_ + a_/c+ 


{a?_ — a\){a\k-hj^ + a^k+h- ^ ^{a^k- + alfe+) 


a+a-^/P 


P 


Ui = 


a?_kj^ — a\k- 
o+/c_ + a_/c+ 


a_/i_ — a+/i+ 7 

a+«- Vp 


(2.25) 


If we substitute these expressions in formulas (|2.23l) . the expressions for the function 
u{x;p) are simplihed and become 


h(x;p) = C±(p)e^^^/“±± 
Here, 


1 


^dioo 


C+{p) = 


Ai/i_ ^ A_/i-i 


2a±^ Jo 
k-j 


[/(0+5(e;p)]e"^'""«'/“±de, ±a;>0. (2.26) 


+ 


‘Ja+y/P 0 ,-y/p Xo{p - a) 


^ t ^ Aih+ , A+/i_ , /c+7 

C-{p) = -7^-^ H-^ + 


2a_^ a+V^ Ao(p + a)' 


A = J_ 1 A = 

^ Ao Vo+ a_ ^ 


k± 


A 

G-j-Aq G+ G_ 


(2.27) 


To finalize our derivations, we apply the inverse Laplace transform and take into consid¬ 
eration the formulas 


Y rc+ioo ^-^a+pt^p ^ 


27ri Jr—i 


^JP 


^/tH 


I {it) 
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Rea > 0, Rec > 0. 


i-c+ioo ^-^a+pt^p 


27 rz Jc—i 


= Erfc 


a 


(2.28) 


Here, Erfc(-) is the the complementary error function. This implies the following repre¬ 
sentations for the function u{x, t) when x is negative: 


u{x,t) = Erfc I — 


+ 


Aoa+ 

1 

2a_\/7rt J- I 


+ 


A, 


^-(a:/a_-^/a+)2/(4t) 


2 \/ta_ / a^y/ni 


fim 




fim 


+ 


a+y/n 


/o Jo 


,-(xla. -5/a+)2/[4h-r)] 9{^:^)d^dT 

y/t^ 


+ 


1 


r r \j^g-(^+g)^/f4a^h--r)] I „-(x-0'^/\Aa'^_it-T)] 9{^ij)d^dT (^2.29) 
Jo J-oo L -I y/t — T 


‘2,0,—^pK JQ J— oo 

For X positive we have 

u{x,t) = --—^Erfc 


+ 


AqU- 

1 

2 a_|_\/ 7rt . 
A_ 

+- 


2y/ia^ 


+ 


A_ 


l-y/ld, J-c 


^-{x/a+-^/a-)^/{4t) 


fim 




1 


2a+v^ Jo Jo 


CL—yTT JO J—oo 

oo r 


Aie-(^+«)V(4<t) ^ g-(x-C)V(4a2 i)l 


(a;/a+—g/a_)^/[4(t—r)] 9iCj T)d^dT 
y/t — T 

2 9{i, T)d^dT 


/[4a^(i-r)] ^-{x-if /[ 4 a/{t-T)] 


y/T^ 


T 


(2.30) 


The total temperature uo{x,t) given by formula (|2.17p is bounded and has different 
limits as X —)• ±oo. When x is kept finite and t ^ oo, the temperature has a finite limit 
independent of x, 

j-a+k- -|- 7 +a_A:+ 


lim uo{x,t) = 


—Ali < X < Af2- 


(2.31) 


t^oo ' AqO+U- 

Here, Xi and X 2 are any finite positive numbers. Formula (I2.3ip is consistent with the 
result obtained by Deconinck et al (2014). If 7 _ = 7 + (7 = 0) and there is no heat source 
{g{x,t) = 0), then the representations (|2.29l) and (|2.30l) generalize to the discontinuous 
case the classical Poisson formula obtained for an infinite homogeneous rod. 

Notice that it is possible to bypass the RHP and derive the representation (I2.26p for 
the function n(x;p) directly by employing the fundamental functions 

^ =-VPl^-«l/a± ^2.32) 


2 a±Jp 

of the differential operators — p. The functions C±{p) are determined in the 

same manner as before from the two conditions (j2.24p . Their expressions coincide with 
those given by (I2.27p . It is evident that the same approach works for any number of 
discontinuities including the case of a finite discontinuous rod with any physical boundary 
conditions imposed at the ends. In this case the fundamental functions (12.321) need to be 
replaced by the corresponding Green functions of the one-dimensional boundary value 
problems. These Green functions are derived in an elementary fashion. If the number of 
discontinuities is n > 3, then instead of two functions C+{p) and C-{p) we have re — 1 
pairs of unknown functions. They are determined by a system of 2re — 2 linear algebraic 
equations following from the 2 re — 2 conditions at the discontinuity points. 


















































2.3 Abrahams-Wickham system of integral equations 

Abrahams and Wickham (1990) analyzed the system 

roo 


u{x) = A / k{x — t)u{t)dt + /(x), 0 < X < oo, 

where the matrix-kernel is given by 


k{x) 




(2.33) 


(2.34) 


A and a are parameters, u{x) = (rti(x),U 2 (x))''~, and /(x) is a forcing vector-function 
prescribed accordingly. To factorize the matrix coefficient of the RHP associated with 
the system (j2.33l) . they expressed the matrix-factors through the solution of a certain 
auxiliary system of integral equations. In the case (|2.34l) that system admits an exact 
solution. In what follows we derive a closed-form solution by a simple method that 
bypasses not only the auxiliary system of integral equations, but also the matrix Wiener- 
Hopf factorization. First we apply the Fourier integral transform to the system (|2.33l) 
and have the following RHP on the real axis 

G{a)U^ {a) = U~ {a) + F'^ {a), —oo<a<+oo, (2.35) 

where 

^ 1 ( a^ + l-2\ -2Ae*““ \ 

+ -2Ae-*"“ ^ 1 - 2A j ’ ^ Jo 

poo /*0 

U~^{a) = / u(x)e*"*(ix, U~{a) = / U-{x)F°‘^dx, (2.36) 

0 J—oo 

the vector-functions U^{a) = {U^ {a){a))^ are analytic in the half-planes = 
± Ima > 0, the vector-function F+(a) = {a),F 2 (a))"'" is analytic in the upper half¬ 

plane and the vector-function u_(x), x < 0, is given by u_(x) = A/q°° k{x — t)u{t)dt. 

Next, instead of factorizing the matrix G{a) we express the function Ui (a) from 
the first equation in ()2.35p and substitute it into the second equation. After obvious 
simplifications this brings us to the new system of functional equations 

a2 + l-2A, , 2Ae*"“, , 

-y—- (a) -" -^1 

+ i + i 

^2 1 _ o\ r)\ —iaa 

+ 1 - + T(o)l + a2 1 _ + -^*<“>1 = ''f <“>' 

Notice that the products (a) and [/]”(«) are analytic in the domains G^ and 

G~, respectively, and vanish exponentially when a —>■ oo and a G G^. For simplicity, 
we analyze further the normal case that is we assume that A ^ [1/4,+oo). Choose 
arg(l — 4A) G (—7r,7r), arg(l — 2A) G (—vr, vr), and denote Aq = \/l — 2A, Ai = y/1 — 4A, 
argAj G (—7r/2,7r/2), j = 0,1. Factorize now the rational functions 

+ 1 - 2A ^ K^{(y) + l-2\ ^ (a) 

g^-|-1 iF|^(g)’ g^ + 1 — 4A Kf(a)’ 
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where 


K^{a) = ( 


a =b 

a zb i J 


T^±t \ (oi±i^o\ 


±1 


(2.39) 


In the general case of the forcing vector-function f{x) we need to introduce new functions, 
'hi (a) and 'I' 2 (a), the Cauchy integrals with the density chosen according to the following 
relations for their limit values on the real axis 


^)''(q:) — (a) = F^{a)KQ (a), —oo < a < -boo, 


^2 («) - ^2 («) = J-) [ + 


F^{a) 2Ae-*““F+(a) 


(a) {a + iXo){a — iXi)' 


—oo < a < -boo. 


(2.40) 


We seek the solution Uj{a) in the class of functions vanishing at infinity. Therefore the 
subsequent application of the generalized Liouville theorem enables us to determine the 
solution of the RHP in the form 


c/r(«) = — 


1 


U^{a) = K^{a) 

1 


'^2 (a) -b 


Koia) [ 
U+{a) = K+ia) 

C2 


Ti(a)+ 

4'+(a) + 


Cl 


OL — iAo_ 

C2 


CX -b zAq _ 
2Ae-““ 


Utia) = 


<(«) 


'hna) -b 


a-b iAo (a-b iAo)(a — *Ai)it'Q (a) 
Cl , 2Ae““i^i+(a) 


Ti(a)+ 


Cl 


ot — iAo 


-b 


a — iAo (a-b i)(a — iAo) 


T+(a) + 


C2 


Ot “h ^Aq 


(2.41) 


Here, Ci and C 2 are arbitrary constants. It is easy to observe that the functions Ui (a) 
and U 2 (a) have inadmissible simple poles at the points a = iAo £ C^ and a = —iAo £ 
C“, respectively. They can be removed if the residues of the functions Ui{a) and 
U 2 (a) at these points vanish. Simple calculations give the following expressions for the 
constants Ci and C 2 : 

c. = (2.42) 


where 


b = 


62 + 1 ’ ^2 ^2 + 1 

2Ae-“^o 


(Ao -b l)(Ao -b Ai) ’ 

di = 2i6Ao4'^(iAo), ^2 = 2i6Ao4')"(—iAo). (2.43) 

The inverse Fourier transformation of formulas (j2.41l) for U^{a) and ?7^(a) yields the 
solution of the original system of integral equations (|2.33n . 


3 Laplace equation in a wedge with mixed boundary con¬ 
ditions: the case of meromorphic functions with periodic 
poles and zeros 

As an illustration of the method in the case when the entries of the matrix coefficient of 
the RHP are almost periodic functions and meromorphic functions with periodic poles 
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and zeros, we consider the following mixed boundary value problem for the Laplace 
operator: 


Att(r, 0) = 0, 0 < r < 00 , 

0 <9 < a, 


^^ 1^=0 = h-{r), 0 < r < ai; 

du 

II 

0 

II 

+ 

2) 

r > oi 

rd9 

u|6)=a =/ 2 -(r), 0 < r < 02 ; 

du 

rd9 

= f 2 +{r), 

9=a 

r > 02 - 


(3.1) 

The function u is sought in the class of functions bounded at r = 0 and vanishing at 
r = oo as \u{r,9)\ < Ar~^, 0 < 9 < a, A = const, 0 < /3 < e, and e is a small positive 
number. This problem can be interpreted as a model of stationary heat conduction if 
(r, 9) = u{r, 9) + T^o is the temperature and Too is the temperature at infinity. Too 
is constant for all 9 E [0, a] and has to be determined a posteriori. The temperature 
?/(*) = Ti is prescribed in the segment 0 < r < ai, 0 = 0 as Ti(r) = /i-(r) + Too and 
in the segment O<r<a2,0 = aas = T 2 (r) = / 2 -(r) + Too. In the rest of the 
boundary, the heat flux q is given: q = /c_/i+(r), r > ai, 9 = 0, and q = /c+/ 2 +(r), 
r > 02 , 9 = a. {k- and are the thermal conductivities of the lower and upper 
boundaries, respectively). It is also assumed that the functions fj+{r) decay at infinity 
as fj+{r) = 0{r-^-^), r ^ oo, j = I, 2. 

3.1 Vector RHP 

Before proceeding with the solution, we extend all the boundary conditions for the whole 
semi-axis as 


u\e=o = fi-{r) + ‘Pi+{r), 


e=o 


= ipi-{r) + fi+{r), 0 < r < oo. 


u\9=a = / 2 -(r) -t <^ 2 +(r). 


du 

rd9 


9=a 


= ip2-{r) + f2+{r), 0 < r < oo. 


(3.2) 


where 


supp/j-(r) C [ 0 ,aj], supp/j+(r) C [0^,00), 

supp(/?j_(r) C [ 0 ,aj], supp(^j+(r) C [0^,00), ( 3 . 3 ) 

and the functions <Pj±{r) need to be recovered. We shall now apply the Mellin transfor¬ 
mation to the Laplace equation and the extended boundary conditions (j3.2l) . Let 

rl poo 

F'fis) = fj{ajp)p^-'^dp, F^{s) = ajJ^ fj{ajp)p^dp, 

pi poo 

^~{s) = Oj Pj-{ajp)p‘'dp, l>+(s) = y^ pj+{ajp)p''~^dp, 


and 


(3.4) 

(3.5) 


poo 

Us{9) = / u{r,9)F~^dr, s E T : Re s = cr E (0,/3). 

Jo 

The function Us{9) solves the differential equation Us{9) -|- s'^Us{9) = 0 and has the form 

Us{9) = Ai(s) cos s9 + A 2 {s) sin s9, (3.6) 

where Ai and A 2 are fixed from the first two boundary conditions in 


as 


= a?[cl>+(s) + Tr(s)], ^2 = 


(3.7) 
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The third and fourth boundary conditions in (13.‘2p constitute the vector RHP 


$+(s) = --G{s)[^-{s) + F+{s)] - F-{s), s E L, 


(3.8) 


where 


Gis) = 


cot as A* CSC as 
\~^ CSC as cot as 


A = 


ai' 


(3.9) 


Without loss of generality we assume that A > 1. The vectors <^=*=( 5 ) = (<h^(s), <l>^(s))''~ 
are analytic in the half-planes : ±Res < ±cj, a E (0,/3), and 0 E . 

3.2 Factorization of the matrix G{s) 

We aim to find two matrices, X'*“(s) and X“(s), analytic in the domains and D~{s), 
respectively, having a finite order at inhnity and solving the following matrix equation 

X+(s) = G(s)X-(s), sEL. 

Denote 

^12(s) \ 

* ’ ■ ( T(o) .T(o )} ' 

Regrouping terms in the same fashion as in section 2.1 we obtain 

A"X2i(s) 


(3.10) 

(3.11) 


Xi, ('S) = -tanasxi (s) 


cos as 


X2j{s) = t?masxtj{s) - 




cos as 


(3.12) 


It is seen that the functions Xi^('S) and X 2 j('®) have simple poles at the points s = 
— ^(n -b E and s = ^(n -|- |) E D~ (n = 0,1,...), respectively, and they have to 
be removed. We seek the functions Xij('S) and xtji^) ™ ^he form of the hyper geometric 
series 

^ V r(as/7r + fc + Z^i)r(fc + t/2) . .ky-kTr/a+ai 

'' r(as/ 7 r + A: + + l/2)r(A; + ^2) ^ ^ 

TW = (3.13) 

^r(-as/7r-bA:-b/3i-bl/2)r(A:-b/32) 

The parameters Vj, ^j, aj, /3j, aj {j = 1,2), and k are to be determined from the 
following conditions: 

(i) the functions Xi 7 ('S) have removable singularities at the zeros of cos as lying in 
the half-plane D~^, 

- sinasxrj('S) + A^X^j(s) = 0, = + n = 0,l,..., (3.14) 

(ii) the functions X 2 ji^) have removable singularities at the points s = ^{n + ^) E D~ 

sin asxtjis) - >G^Xij{s) = 0, ® ^ ’ n = 0,1,..., (3.15) 
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(iii) the functions Xij('S) niay have simple poles at the zeros of tanas lying in D"*", 
s = —vrn/a, n = 0,1,...; it follows from (13.121) that then the functions XyC'S) have 
removable singularities at these points, 

(iv) the functions xtji^) have simple poles at the points s = Trn/a G D~, 
n = 1, 2,.... Then the functions X 2 ji^) have removable singularities at these points. 

The conditions (i) and (ii) when satisfied give 


/?2 — 1 , /^2 — 1 , / 3 l — 1 — fJ - l , 02 — 1/2 + 1^1 — // i , 

1/2 = ai +/xi - 1 / 2 , K = a 2 = ai + -ll. ( 3 .I 6 ) 

2a 

The other two conditions, (iii) and (iv), determine i^i and ai: z/i = 0 and ai = 1. 
Without loss of generality, ai = 0. This brings us to the following one-parametric family 
of solutions: 

_ . . _ ^ r(ag/7r + k)T{k + ni + l/2)(-l)^A-^^/» 

^ ^ r(as/7r + A;-I-+ 1/2)A;! 




(^) = E 

k=0 


Ti-as/n + k + l)T{k - m + l/2)(-l)^-^i 
r(— as/vr + k — fii + 3/2)kl 


(3.17) 


On choosing /ii = 0 for j = 1 we obtain the functions Xu ('5) xJl(' 5) vanishing at 

infinity as xri('®) ~ 0(s“^/^), s G D~, and yore's) = 0(s“^/^), s G D~^, and admitting 
the following representations through the beta and Gauss functions: 


_ /as 1\ /as 1 as 1 

XnW = B(^V2j'^lV'2’V+2--^ )■ 


as . , 1\ ^ / as „ 1 as , 3 .-jr/o 


X^,(S) = ^ K + i-A 

TT Z/vTT Z TT Z 


(3.18) 


The function xri(®) analytic in the half-plane D~ and meromorphic in , while the 
function xJl(®) analytic in and meromorphic in D~. In the case j = 2, we put 
= 1 and derive the other two functions 


X22('S) = 


- / A „ /as 3\ ^ /as 3 as 3 ,_^/„A 

W») = b(-,2)^(v2-V+2--" ')• 

G “ + 1 , _ T f G “ + 1 , _ 1 52 + 1 _a-.a 

V7r Z J \ TT Z TT Z 


(3.19) 


Notice that the function Xio(s) vanishes at infinity, while the function Xoo(s) grows as 

s ^ 00: xr2('S) = 0(s“3/2)^ ^ <9^51/2^^ 

The functions Xi/ (s) and X 2 j(^) expressed through the functions ()3.18p and (I3.19P 
by formulas p3.12p . They are analytic in D'^ and D~, and their order at infinity is 
determined by the order of Xi/ (s) and x^ (s), respectively. This completes the exact 
factorization of the matrix G(s). The factorization we found is not unique: by choosing 
the parameter /j,i as an integer different from 0 and 1 we can construct another set of 
the functions xr 2 ('®) X 22 i^) that is why another factorization different from the 

one found. The Wiener-Hopf matrix factors X+Ps) and X~ given by (13.181) and (I3.19P 
have the following asymptotics at infinity: 


X±(s) 


0(s-V2) 0(s-3/2) A 

0(S-V2) 0(,l/2) j’ 


(3.20) 
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3.3 Exact solution of the heat conduction problem in a wedge 

On having factorized the matrix G{s) we replace G{s) by the product 
in the boundary condition (13. 8p of the vector RHP 

[X+(s)]-i[<l>+(.) + F-{s)] = -i[X-(.)]-i[cl>-(.) + F+{s)], s E L. (3.21) 

Represent next the functions fi-{r) and / 2 -(r) as 

fj_{r) = f,+T*{r), fj = -T^ + TJ, j = l,2, (3.22) 


where T? +Tj (r) = Tj{r) is the prescribed temperature in the lower boundary 0 < r < ai 
and the upper boundary 0 < r < 02 of the wedge, T° = const, and (r) = 0(r'^^}, r —)• 0, 
jj > 0. We also introduce the Mellin transforms 

f-(s) = £ T;(ajp)p^-^dp, j = 1, 2, (3.23) 

and the Cauchy integrals 




1 f [X+{(7)]-^f-{a)da 
2-Ki Jl a — s 




1 f [X (u)] ^F'^{a)da 


(3.24) 


a{a — s) 

Employ now the continuity principle and the Liouville theorem and derive the following 
representations of the solution: 


$+(s) = (-/ + X+(s)[X+(0)]-^) ^ - X+{s)[^+{s) + n+is)], 


$-(s) = -X-{s)[X+iO)]-^f + sX-{s)[^-{s) + n-{s)], (3.25) 


where T = {Ti^T 2 )^■ Further, since the functions x^(s) are growing at infinity, and the 
matrix factors have the asymptotics (13.201) . the second components of the vectors $”( 5 ) 
and <h+(s), in general, have the asymptotics <^ 2 ( 5 ) = 0(s^/^), $^( 5 ) = 0(s“^/^). Due 
to the Tauberian theorems for the Mellin transforms this causes an infinite temperature 
as r —)• at the lower {j = 1) and the upper {j = 2) boundaries. Also, as r —>■ aj, 

the heat flux has a nonintegrable singularity of order 3/2. To derive the condition which 
guaranties the boundedness of the temperature and integrability of the heat flux, we 
denote the limits 




(D?,D°2)^ 


lim s'h^(s), 

s^oo,sGD^ 

lim sD^(s), 

s^cxd,sGZ)=*= 


(3.26) 


and compute [Al+(0)] ^T. It is directly verified that 


Xn(0) = -vr + xJl(O), X^ 2 ( 0 ) = -tt + X^ 2 ( 0 )> (3-27) 


and 

X^i(O) = 2tan-i ^+^( 0 ) = + 2tan-^ ^3 28) 
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Continuing, the vectors <h^(s) have the asymptotics at infinity we need if and only if 


- Ti + ^2 



A+(0) 

X^i(O) 


+ ^ 2)1 


(3.29) 


where A+(0) = detX+(0) = —If this condition is fulfilled, then as s — 00 , 
<I>“(s) = s E D~, and <I)+(s) = 0(s“^), s E D~^. The condition can be easily 

satisfied by the corresponding choice of the parameter Too undetermined at this stage. 
Since Tj = —Too + ^/) we transform equation (|3.29p and find the temperature at infinity 

Too = - tan"^ (r° - T^) + + 0^). (3.30) 

TT 

In particular, if T*{r) = 0 and fj+ir) = 0, then = 0 and = 0- assuming 
further that T^ = 0 (the case considered by Abrahams and Wickham (1990)), we obtain 
the simple explicit formula 


--tan-^ ^ 1 


T2°. 


(3.31) 


Finally, if oi = 02 , then A = 1 and therefore Too = 5^2 • The same result for the case 
A = 1 was derived by Abrahams and Wickham (1990). 

We wish also to recover the function u{r, 0) and study its behavior as r —>■ 00 and 
r —)• 0. On applying the inverse Mellin transform to ()3.6I) and using (13.71) and (13.Sp it is 
possible to have 


M(r, 0) 



[<1>;^ (s) + T]+(s)] cos(a 


e)s 



— S 


+ [^2 ('S) + ^2 (s)] COS Os — 

\a2 


ds 


s sm as 


(3.32) 


Denote by kj those singular points of the functions (s) in the domain D which have 
the smallest real part among the singular points of F^{s) in D~. Then we can derive 


ti(r, 9) = 0(r ^), r — 00 , /3 = min{ 7 r/Q:, Ki, K 2 } > 0. 

In particular, when /i+(r) = / 2 +(r) = 0 and r —>■ 00 , 


u[r, 9) = — cos — 
vr a 




—Trjcx. 


I -) ( — 

aj V 02 


— IT jo. 


+ 0{r 


—2tt/ cx.\ 


(3.33) 


(3.34) 


and as 0 = a/2, u ~ r 00 , c is a nonzero constant. 

If we want to determine the behavior of the function u as r —>■ 0, we need regroup 
the terms in the integrand (13.321) in order to get rid of the functions d>^(s) and 4>2 (s). 
We have 

uir, ^ I (~) 

+ [cI>+(.) + T2-(s)]sin0s(^) 1^. (3.35) 
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Accepting the representations (I3.‘22p by the Cauchy theorem we derive 


u{r, 6) 


1 - -) + -T^ - Too, 

a J a 


0 . 


(3.36) 


Finally, we notice that the solution we obtained have the following asymptotics at the 
points where the type of the boundary conditions is changed: 


u ^ Cj, r ^ a- 


de ^ 


e = 


0 = 0, 3 = 1, 
0 = a, j = 2. 


(3.37) 


4 Helmholtz equation in a strip: nonperiodic poles and 
zeros 


With the example of the Helmholtz equation in a strip with a cut inside we shall show how 
the method of the RHP can be generalized when the zeros and poles of the meromorphic 
entries of the RHP matrix coefficient are not periodic. Consider the Neumann boundary 
value problem for the Helmholtz equation in the doubly connected domain H \ S' with 
H being a strip, H = {|xi| < oo, —bt < X 2 < &+}, and S being a cut, S = {0 < xi < 
a,X2 = 0 =^}, 

(A + A:o)«° = 0, (xi,X2)en\S, 

-— = 0, lxi|<oo, X 2 = 

0X2 

du° 

-— = /°(xi), 0 < xi < a, X 2 = 0=^. (4.1) 

0X2 

This problem can be interpreted as an antiplane problem for a strip H with a crack 
S when the strip boundary is free of traction and the crack faces are subjected to the 
oscillating load r 23 = , t is time, uj is the frequency, and G is the shear modulus. 

This problem also describes acoustic wave propagation in a guideline with an acoustically 
hard screen S inside. 

We assume that Imfco >0. R will be convenient to work in the dimensionless 
coordinates x = xi/a, y = X 2 /a. Denote k = ako, b± = b’^/a, f{x) = af°{ax), u{x,y) = 
u°{ax,ay). This transforms the problem to 


(A + k‘^)u = 0, (x, y) G {|x| < oo, —6_ <y < 6+} \ {0 < x < 1, y = 0^}, 


du 

dy 


= 0, |xl < oo, y = ±6±, 


du 

— =/(x), 0 < X < 1, y = 0^, ti(x, 0^) — n(x, 0“) = (/7i(x), |x| < oo, (4.2) 

where suppy^i C [0,1] and (pi{x) is an unknown function in the segment [0,1]. 

On applying the Fourier transform with respect to x 


Ua{y) = / u(x, y)e*"'^dx 

J — OO 

to the Helmholtz equation we easily obtain 


(4.3) 


Ua{y) 


Gf cosh 7 (y — 6_|_) + sinh 7 (y — 6+), 0 < y < 6+, 
Cf cosh 7 (y + 6_) + sinh 7 (y + 6_), —6_ < y < 0, 


(4.4) 
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where 7 = y/— k'^ is the single branch of the algebraic function 7 ^ = fixed by 

the condition Re 7 > 0 or, equivalently, 7 ( 0 ) = —i in the a-plane cut along the straight 
line joining the branch points k and —k and passing through the infinite point. The 
boundary conditions on the sides y = and the discontinuity of u and continuity of 
its normal derivative in the line y = 0 yield 


+ _ _ , _ sinh 76 -^j^(a) _ sinh{a) 

^ ^ ^ sinh 7 ( 6 ++ 5_) ’ ^ sinh 7 ( 6 ++ 6 _) ’ 


(4.5) 


where 


$+(«) = / 
40 


(4.6) 


To derive the RHP, we hrst extend the Neumann boundary condition on the line y = 0 


du 

dy 


f-{x) + y:>2-{x) + ip2+{x), -00 < X < 00 , y = 0 =^, 


(4.7) 


where f-{x) = f{x) if 0 < x < 1 and it vanishes otherwise, <p±{x) are unknown functions 
such that supp(?!) 2 - C (—oo,0] and supp 02 + C [l,oo). In order to apply the Fourier 
transform to the boundary condition ()4.7I) . we introduce the following integrals 


^■*■( 0 :) = [ /(x)e*“*(ix, 

Jo 

/ O POO 

<y 92 _(x)e*"®dx, <i>^(a)= / ip 2 +{x + l)e"^'^dx. (4.8) 

-00 J 0 


We assert that <hj’'(a) and F~^{a) are entire functions which are analytic in the upper 
half-plane and have an essential singularity at the infinite point in the lower half-plane 


and also 

4>+(a) = e*“4.r(a), 

F+{a) = F^F-{a), 

(4.9) 

where 

rO 

(a) = J <y9i(x -b l)e*“^(ix. 

F~{a) = J f{x + l)F^^dx. 

(4.10) 


This implies that the boundary condition (j4.7jl is equivalent to the following vector RHP: 


‘h^(a) = G(Q:)<h (a) -|- F (a) J, —00 < a < -|-oo. 


(4.11) 


Here, 



g{a) = ysinhyb^ sinh 76 _csch 7 ( 6 + -|- 5_). 


(4.12) 


The new feature here is that the zeros and poles of the meromorphic entry of the 
matrix G{a), the function g{a), are not periodic, and we will pursue the RHP in a 
way different from the one proposed in the previous section. First we factorize the 
meromorphic function g{a) 

9(a) = (4.13) 

K [a) 

where 

K^{a) = {a± k)'^^^'^K^{a), —00 < a < -|-oo. 
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Ko(a) = exp — / 

, vrz Jo 


a 


log 


sinh 76 + sinh 76 - d(3 


a ^(— 00 ,+ 00 ), (4.14) 


sinh7(6_|_ + 6_) (3“^ — 

and K^{a), —00 < a < + 00 , are defined by the Sokhotsky-Plemelj formulas. Then we 
rewrite vector equation (14.IIP as 




<^2 (a) T + / \ 4)i (a) 

^ ^ Tlfal =- l-^-L - 

K+{a) 


e *“<h 2 (a) 


K (a) K {a)g{a) 




where 


4'_(a) = 


K-{P)F+{P)dp 

(3 — a ’ 


T+(a) = 


.00 F-{/3)d/3 

Loo K+{(3){P-a)- 


(4.15) 


(4.16) 


The right-hand sides of the first and second equations in (|4.15l) are analytic everywhere 
in the domains C~^ : Ima > 0 and C~ : Ima < 0, respectively, except at the zeros of 
the function < 7 ( 0 ). Denote by am± the zeros of the functions sinhy^^ lying in the upper 
half-plane, 


Q^mdi 


= Jk^- 




Im am± > 0 , m = 1 , 2 ,.... 


(4.17) 


Then in C'"*', the function g{a) has simple zeros at the points a = k and a = am±, 
m = 1,2,.... In the lower half-plane, the zeros are a = —k and a = —am±-, m = 1, 2,.... 
In order to remove the simple poles of the right-hand sides of equations (j4.I5jl . we 
introduce the functions 


4+ °° 


D (a) = 


a + k 
^0 


/ 1 + 4 + 

^m-\- _j_ -^m- 


~=1 \ ^ ^m+ ^ H“ ^m— 


+ E 


A, 


•m+ 


-7 




Q; k rn=l ^m+ ® Q^m— 


(4.18) 


Assume that the residues of the functions D“(q;) and D'’'(q;) at their poles are chosen such 
that, when they are subtracted from the left and right-hand sides of the first and second 
equations in (14.151) . respectively, the poles are removed. On employing the Liouville 
theorem we find the solution of the RHP 

^^(«) = K-{a) ’ ^ K+{a)[-^X{a) + H+(a)], 


^i{a) = K (a)[^+(a) - H+(a)] -— 


K {a)g{a 


■[4'l(a)-hH (a)], 


, Tl(a) + H-(a) e*“i7+(a)._, , 


K+{a) 


aioi) 


(4.19) 


It is directly verified that ‘h)^(a) = (a). The coefficients Aq , and A*_ need 

to be fixed by the conditions 


res <hi*^(a) = 0, res <ht(a) = 0, res (a) = 0, m = l,2,.... (4.20) 

Q;=i/L —^mziz 0 !.=— 
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These conditions if satisfied guarantee that the zeros of the function g{a) are removable 
singularities of the functions <l>^(a). On calculating the residues we rewrite the conditions 
(I4.2U|1 as a system of linear algebraic equations 




/i± °° / \ 

2k m=l \ ^ ^m+ k T CXm— j 


±5oe^^^i{±k), 


A-. + 


A+ 

° +E 


4 + 

-^m+ 


+ 


4 + 
^m— 


Oin± T k Yn=l \ Q^n± T CX-m— 


= 6n±e^^-^A>+{an±), 




*an± 


An 


+ E 


A: 


m+ 


A: 


CUnzt T k \ ®n± T 0'm+ ®n± T Oim— 

=-5„±e*“"±TI(-an±), n = l,2,.... 


(4.21) 


Here, 


Sn± = 


(-1)"-(1 + k/an±)[K^{an±)f sm[7m{b+ + b-)/b±] 


b± sin(7rn5=p/6±) 


5o = 


b, ^ 6 _ 


[K+{k)f. 


(4.22) 


If we assume the symmetry of the problem, 6_ = 6+ = 6, then the formulas can be 
simplified. Indeed, g{a) now becomes ^ytanhy^, and therefore its zeros are given by 
am = y/k'^ — Tr'^m^/b'^ (oq = k), Imam > 0, m = 0,1,.... The functions H^(a) are 
represented by the series 

Si±(a) = f; (4,23) 

^ a± am, 


m=0 


and the infinite systems reduce to 


m=0 


AT 


= ht, n = 0,1,.... 


(4.24) 


Here, 

ht = T<5ne*“"'&^(T«n), 5n = [K^{a^t. (4.25) 

a^iO 

Notice that in both cases, 6_ ^ 6+ and 6_ = 6+, the unknowns A+ and A~ exponentially 
decay as n ^ oo, |AJ| < c± and q are some positive constants. 


5 Conclusion 

We have developed further the algorithm for the vector RHP when its matrix coeffi¬ 
cient entries are meromorphic and almost periodic functions. In the simplest case, when 
the meromorphic functions have a finite number of poles and zeros (rational functions), 
regardless of the dimension of the vector RHP, the exact solution can always been con¬ 
structed. We have illustrated this approach by solving the inhomogeneous heat equation 
with a piece-wise constant diffusivity. In the case of two discontinuities, 0 and oo, we 
have derived a simple representation of the solution that generalizes for the discontinuous 
case the classical Poisson formula for an inhnite rod. 

We have also shown that, when the meromorphic functions have periodic zeros and 
poles only, then it is possible to derive the Wiener-Hopf factors of the RHP matrix 
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coefficient in a closed form in terms of the hyper geometric functions. This technique has 
been employed for finding the exact solution of the vector RHP associated with a mixed 
boundary value problem for the Laplace equation in a wedge when on finite segments 
of the wedge sides, the function is prescribed, and on two semi-infinite segments, the 
normal derivative is known, and no symmetry that would allow for decoupling of the 
problem is assumed. 

In the general case, when the zeros or poles are not periodic, the solution can be 
derived by quadratures and some exponentially convergent series. The series coefficients 
solve an infinite system of linear algebraic equations of the second kind whose rate of 
converge is exponential. This technique has been used for solving the Neumann boundary 
value problem for a strip with a finite slit inside parallel to the strip boundaries. 
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